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The stationary state solutions of the Ising-type thin films with difi^erent layers in the presence of 
an external oscillatory field have been examined within the effective-field theory. The exhibition 
focuses on understanding of the external field frequency and amplitude effect on the overall behavior 
of system for several powerful treatments. The particular attention has been paid on the dynamical 
parametric evolution of the special point at which the critical temperature of the three-dimensional 
infinite bulk system, where the surface and the modified exchange parameter are of no importance. 
Important characteristics of the films, such as surface enhancement phenomenon, effect of thickness 
on the dynamic process etc. have been introduced together with well known other characteristics. 
An attempt has been made to explain the relations between the competing time scales (intrinsic 
microscopic relaxation time of the system and the time period of the external oscillatory field) and 
frequency dispersion of the critical temperature coordinate of the special point. 
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I. INTRODUCTION 

The conceptual frame of surface enhancement phe- 
nomenon on the finite magnetic materials, especially on 
semi-finite systems and thin films, has not lost its cur- 
rency yet and recently there has been growing interest 
for this mature subject of solid state physics because a 
variety of apparently very different phenomena can be de- 
scribed by both experimentalists and theorists. Thin film 
of nanometer thickness has been fabricated by various 
coating techniques and has many applications in tech- 
nology. Advances in experimental techniques today allow 
the crosschecking of such far-reaching predictions. In the 
past quarter century, some preliminary experimental ev- 
idence has been presented for the existence of a surface 
type of magnetic ordering since the magnetic properties 
of free surfaces drastically differ from the bulk material, 
because the free surface breaks the translational symme- 
try, i.e. surface atoms are embedded in an environment 
of lower symmetry than that of the inner atoms and con- 
sequently the exchange constants between atoms in the 
surface region may differ from the bulk value. 

One of the early experimental studies involving 
ground-breaking discoveries in the understanding of 
surface-enhanced magnetic order and magnetic surface 
reconstruction on Gd(OOOl) by Weller et al. dating 
back to the last quarter of the century, there is an empir- 
ical inference that the 4/ spins of the surface are not fer- 
romagnetically coupled to the bulk moments. They have 
presented spin-polarized low-energy-electron diffraction 
(SPLEED) and magneto-optical Kerr effect (MOKE) 
studies performed in situ that demonstrate the existence 
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of surface-enhanced magnetic order. Their results re- 
vealed the necessity of new Monte Carlo (MC) calcu- 
lations for systems with antiferromagnetic (AF) perpen- 
dicular surface layer coupling. Also, Diirr et al. ^ have 
shown in the same period that the temperature depen- 
dence of the long-range order parameter in thin (1-3 ML) 
films of bcc Fe on Au( 100)— including one monolayer—, 
as measured with SPLEED and spin-polarized secondary 
electron emission spectroscopy (SPSEE) as a comple- 
ment, represent a physical realization of a truly two- 
dimensional phase transition. Additionally they reported 
that the thickness independence of the critical behavior 
confirms the universality hypothesis. In the former of the 
two unassailable experimental observations published in 
the early of two decades ago by Li and Baberschke [3] , in 
situ magnetic resonance (MR) measurements of ultrathin 
Ni(llO) prepared in ultrahigh vacuum have been carried 
out to determine the thickness dependence of the critical 
exponent /3 and the Curie temperature Tc- They con- 
cluded that the existence of a crossover in the magnetic 
properties of Ni films as a function of thickness and the 
critical exponents determined in their experiment agree 
very well with the known exponents for a 2D Ising-type 
system and bulk Ni. Subsequently, thickness-dependent 
Curie temperature of 5-100 monolayers Gd on W(IOO) 
and its dependence on the growth conditions are deter- 
mined by in situ ac-susceptibility (xac) measurements by 
Farle et al. Q in the second one. They found that the 
Curie temperature of carefully prepared layer- by- layer- 
grown Gd(OOOl) films varies from rc(bulk) = 292. 5K 
to Tc(5ML) — 120K and it was also demonstrated that 
different growth conditions change the film's Curie tem- 
perature dramatically. After the aforementioned seigno- 
rial studies, a comprehensive examination propounded in 
the subsequent years by Poulopoulos and Baberschke 0] . 
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They discussed the fundamental magnetic observables, 
i.e. magnetic moment per atom, Curie temperature, sus- 
ceptibility and magnetic anisotropy, for idealized proto- 
type thin films like Fe, Co, Ni on metal substrates such 
as Cu, W, Re. They also presented in depth satisfactory 
knowledge about studies on trilayers, i.e. magnetic thin 
films separated by a spacer, like Cu in the relevant com- 
prehensive work. Apart from these, important studies 
are so located as well as experimentally, the interesting 
phenomenon of surface enhanced magnetic order, i.e. the 
coexistence of an ordered surface with a disordered bulk 
has been observed in two of the 4/ rare-earth metals: 
Gd and Tb. Former experimental measurements sug- 
gest the possibility of large surface anisotropics in thin 
films [1, H, H, HI]. Beside the relevant elements, surface 
ordered magnetically can exist also for Cr [13, [HI. In 
general, one finds that due to the reduced coordination 
number, the critical temperature is lower at surface in 
thin films and decreases with decreasing film thickness 
[l3 - [l3 |. As presented by Detzel et al. [l5| and Stam- 
panoni et al. [l^, for fee Fe films on Cu(lOO) substrates 
the Curie temperature increases strongly from one to two 
monolayers, but decreases in thicker layers. For the Fe 
films grown on Ni/Si substrates, it has been observed that 
a transition between magnetic and nonmagnetic phases 
occurs at a critical thickness [T7l - l20| . With the recent ad- 
vances of modern vacuum science, in particular epitaxial 
growth techniques (especially modern beam epitaxy) it 
is now possible to grow very thin magnetic films of con- 
trolable thickness and this has stimulated renewed inter- 
est in both experimental and theoretical film magnetism 

MM- 

Up to date theoretical treatments on the semi-infinite 
systems and thin films in the literature have continued 
simultaneously with experiments. Since the early experi- 
mental observations, several theoretical and/or computa- 
tional studies on pure crystalline thin films have been put 
forward to explain nature of surface enhancement mecha- 
nism. First, by Binder and Hohenberg [23i], Ising models 
with modified exchange interaction on the surface are 
considered. The critical value of the modified exchange 
for surface ordering was found from high-temperature- 
series expansion (HTSE) and compared to mean-field 
theory (MFT) value in their pioneering work. They re- 
ported that there is a temperature region in which the 
surface behaves like a bulk two-dimensional Ising model 
near its phase transition for values greater than the criti- 
cal value of relevant exchange. The critical exponents ex- 
perience a crossover at the critical value of the exchange, 
which is reflected in poorly behaved series, and effective 
exponents differing from the true ones for values less than 
the critical value of the exchange. The average magneti- 
zation of the simple cubic Ising films (with different thick- 
nesses atomic layers) is calculated as a function of tem- 
perature using the MC technique by Binder in the same 
year [2^. Computational details have been also given 
by the author. The semi-infinite Ising model with an 
arbitrary number of surface magnetic couplings different 



from the bulk exchange constant, is solved in the MFT by 
Aguilera-Granja and Moran-Lopez in the middle of three 
decades ago [23. Cardinally they showed that the exact 
expressions for the critical couplings leading to a surface 
Curie temperature higher than the bulk one. Landau and 
Binder [2g extended the MC simulations and they pre- 
sented the results of phase transitions and critical behav- 
ior at the surface of a simple cubic Ising model. The work 
contains the detailed extraction of surface and bulk prop- 
erties as a function of temperature and surface coupling. 
Also, the surface-bulk multicritical point is located with 
improved precision and crossover behavior is studied. In 
the maturation stage of the literature, review of surface 
magnetism was given by Kaneyoshi [27j . In particular the 
interplay of magnetization and anisotropy at a surface is 
discussed, after a survey of experimental and theoretical 
results for magnetic moment and anisotropy at surfaces 
or interfaces of semi-infinite magnets and thin films in the 
article. Ferroelectric films have been studied to obtain 
the polarization and the dependence of the Curie temper- 
ature on the long-range exponent in two works [2^,[2^. A 
cubic ferroelectric lattice is assumed to consist of pseudo- 
spins with interactions, parallel study using a distorted 
lattices, such as tetragonal lattice, gives the same depen- 
dence of the Curie temperature and polarization accord- 
ing to Wang et al. [28]. Yao et al. [2^ presented that 
the long-range interaction makes the transition temper- 
ature higher. In both studies, films have been described 
by transverse Ising model (TIM) with long-range interac- 
tions. The phase transitions of in a transverse spin-1/2 
Ising film has been also examined in detail within the 
framework of the effective-field theory (EFT) by Saber et 
al. 30] . It has been found that if the ratio of the surface 
interaction to the bulk one less than a critical value, the 
critical temperature of the film is smaller than the bulk 
critical temperature and as the film thickness is increased 
further, critical temperature increases and approaches 
asymptotically to the bulk-critical temperature for large 
values of the thickness. Additionally they reported that 
the critical temperature is larger than both the bulk and 
the surface critical temperatures of the corresponding 
semi-infinite system and as the film thickness is increased 
further, critical temperature decreases and approaches 
asymptotically to the bulk one, for large values of the 
thickness, the surface magnetic transitions observed for 
greater values of the critical value of modified exchange 
in the corresponding semi-infinite system. Also, the layer 
magnetizations and their profiles have been presented as 
an essential detail and they illustrate the existence of one 
defined critical temperature of the film. In recent years, 
Laosiritaworn and co-workers [3l| have used MC simu- 
lations and the MFT to observe the magnetic behavior 
of Ising films with cubic lattice structure as a function 
of temperature with thickness and they found that the 
magnetic behavior changes from two-dimensional to the 
three-dimensional character with increasing film thick- 
ness and this finding made great strides in semi-infinite 
concept. Both the crossover of the critical temperature 
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from a two-dimensional to a bulk value and shift expo- 
nent are observed. The same computational technique 
has been used by Cossio et al. [32] to explain the tem- 
perature dependence of the magnetization, the magnetic 
susceptibility and also the fourth-order Binder's cumu- 
lant. Ferromagnetic (F) Ising-type thin film in the pres- 
ence of a time dependent external oscillatory magnetic 
field has been investigated to model the hysteretic behav- 
ior of the system by Laosiritaworn [s^l • In the work, MC 
simulation was used and thickness dependence of hys- 
teresis properties with varying frequency and amplitude 
of the external field has been presented by the author. In 
MC simulations on semi-infinite systems, a linear dimen- 
sion parameter should be much larger than the other two 
to maintain the quasi two dimensional structure of the 
films, so it has been taken in consideration to preserve the 
films geometry. The power law scaling relations among 
the hysteresis properties, the thickness and the field pa- 
rameters have been found i.e. scaling of the hysteresis 
loop (HL) area performed and scaling exponents were 
also reported. However, the rich behavior— which was 
presented in literature casting above— depending on the 
dynamic parameters has not been covered in the article. 
In light of those results, it appears that hysteresis prop- 
erties of thin films under external perturbation can be 
predicted which provide another successful step in mod- 
eling ferromagnetic materials. In a brand new paper, 
the effect of the random magnetic fields distributed by a 
Gaussian distribution centered at zero on the phase di- 
agrams and ground state magnetization of the thin film 
described by TIM has been investigated by Akmci [33 |. 
Particular attention has been paid on the evolution of 
the special point coordinate with distribution parameter 
in relevant study. The author found that rising the dis- 
tribution parameter makes no significant change special 
point in both planes, but gives rise to a decline in co- 
ordinates of it and after a certain value of distribution 
parameter, the special point disappears in both planes 
simultaneously. 

In the semi-infinite systems, depending on the ratio 
between surface and bulk exchange interaction, the sys- 
tem may order on the surface before it orders in the bulk 
which is called extraordinary transition. In the contrary 
of this, ordinary transition means that the surface critical 
temperature is the same as the bulk transition temper- 
ature. There is an exact consensus in literature that, 
the intersection point between these two transitions is 
called 'special point' (in other words 'crossover point'). 
As the film gets thicker, it approaches the semi-infinite 
system, and this unusual effect shows itself as an intersec- 
tion point in the phase diagrams plotted in critical tem- 
perature versus surface exchange interaction plane for the 
films that have different thicknesses. 

It is well known that the EFT is one of the most pow- 
erful methods that determines the meniscus which sepa- 
rates several phases in the relevant planes, based on the 
use of rigorous correlation identities as a starting point 
and utilizes the differential operator technique firstly de- 



veloped by Honmura and Kaneyoshi [sg!]. Although the 
conventional version of the method fails to find an ex- 
pression for the free energy, since it takes into account 
the self spin correlations, the method is superior to MFT 
which neglects the thermal fluctuations via neglecting the 
self spin correlations. Thus, it is expected from EFT to 
obtain more reasonable results than MFT for these sys- 
tems, as in the case of static Ising model. Therefore in 
this work, we intend to probe the phenomena mentioned 
above of Ising-type semi-infinite systems and thin films 
in the presence of a time-dependent external oscillatory 
magnetic field by using the EFT. These types of pertur- 
bations constitute an important role in material science, 
since the time-dependent external effects may reveal the 
origin of some important macroscopic behavior pattern 
which is still open for inspection in intrinsic surface mag- 
netism. Eventually, for these purposes the outline of the 
article is as follows: We briefly describe the formalism 
and the method used for the system in Sec. 2. Numeri- 
cal results and discussions are summarized in Sec. 3, and 
finally Sec. 4 contains remarks about our conclusions. 



II. METHODOLOGY 

To investigate the dynamical transitions, surface en- 
hancement phenomenon, dynamical symmetry loss (sim- 
ilar to the phrase 'symmetry breaking' in the static 
model) , strength of the magnetic order in different layers 
and many other dynamical features, one simple choice 
may be a dynamical Ising-type simple cubic isometry 
which has an inner coordination number 2 = 4 defined 
on the three dimension with a time dependent external 
oscillatory (in time but uniform over the space) magnetic 
field studied by EFT. For this purpose, we consider the 
following Hamiltonian, 

H = ~^J^JS^Sj-h{t)^S^ (1) 

where Si is the spin operator at a lattice site i and any 
spin variable can take the values Si = ±1. As is known, 
(...) subscript bracket symbolizes the nearest neighbor- 
ing in first summation. The second summation is over all 
the lattice sites. The exchange interaction between 
the spins on the sites i and j takes the values according to 
the positions of the nearest neighbor spins. Two surfaces 
of the film have the intralayer coupling Ji . The interlayer 
coupling between the surface and its adjacent layer (i.e. 
layers 1, 2 and L, L — 1) is denoted by J2. For the rest 
of the layers, the interlayer and the intralayer couplings 
are assumed as J3. The system has three exchange in- 
teractions where Ji, J2, J3 > favors a ferromagnetic 
leaning of the adjacent sites as shown in Fig. ([T]) and the 
Zeeman term describes interaction of the spins with the 
field of the sinusoidal form 

h{t) ^ hocos{ujt), (2) 
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where t is the time and ho is the amphtude of the oscil- 
latory magnetic field with an angular frequency w. 

Conventionally, the layer agents remaining invariant 
under the appropriate symmetry operation in a bulk sys- 
tem are defined with certain criterions: Spin values, types 
and number of it's neighboring interactions. The relevant 
agents depict thermo-magnetical behavior of whole ma- 
terial whereas in a little more realistic systems (such as 
magnetic thin films) , one should construct the equations 
of state for each layer agents to control the dependency 
of the characteristics with the film thickness hence, such 
criterion-related validity is lost. Our system is in contact 
with an isothermal heat bath at given temperature T. 
So, the dynamical evolution of the system may be given 
by non-equilibrium Glauber dynamics [11] based on a 
master equation. The dynamical equations of motion for 
each layer are in the form of 

rj^{4'^) = -{4'^) + (tanh[/?(£;f ) + hit))]), 
rj^{sf^) = -{sf^) + (tanh[/3(£;P + h{t))]), 

rj^{4'^) = -(4'^) + (tanh[/3(£;f-) + hit))]) 

rj^{4'''^) = -{4'^'^) + (tanh[/3(£;f + hit))]) 
^^(4"^^) = -(^i^^) + (tanh[/3(i?P + hit))]), (3) 



.V 




FIG. 1. (Color online) Magnetic unit cell for the simple cubic 
thin films of two difi'erent thicknesses as L = 4, 5. The lattice 
has a finite size only in z direction as other two directions are 
defined on the interval such that x,y e (— oo, oo). Green, ma- 
genta and blue cylindrical edge rendering functions represent 
Ji , J2 and J3 ferromagnetic exchanges respectively, as well as 
spherical vertices symbolize the magnetic items. 



with 

5=1 

p(2) _ 7 Jl) , 7 J2) J (3) 

5=1 



5=1 



p(i-l) _ 7 , Y^ 7 AL-l) J (L) 

A - J3SS' +2^J3Sg +J3Ss> , 

5=1 

=^2.r^^+E^i4"^ (4) 
5=1 

where l/r is the transition per unit time in a Glauber 
type stochastic process, (3 — l/fc^T and fcs represents 
the Boltzmann constant. The subscript and superscript 
in a spin variable symbolizes the lattice site and layer 
number respectively as well as £'{'"^ is the local field act- 
ing on the site i in fcth layer. 6 and d' subscripts in 
the energy contributions represents in-layer and out-of- 
layer nearest neighboring and this type of indication also 
is reinforced by superscripting with regard to control the 
relation between production index and the nearest neigh- 
boring argument in progressive aspects. The average 
brackets in Eq. ([3]) stand for the usual canonical thermal 
average. 

In order to handle the second term on the right-hand 
side of the Eg . dp one can use the differential operator 
technique [H, [STf ? By using the technique, Eq. ([3]) gets 
in the form of 

^) = ~{s\'^) + (exp(i?f V))/(x)U=o, 
Tj^is?^) = -{4'^) + (cxp(i?f V))/(x)U=o, 

rj^isl'^) - -isl'^) + (exp(i?f V))/(x)U=o, 

-|(4""'^) = -{s\'-'^) + (exp(£;f -^V))/(x)|.=o, 

t|#^) = -(^f') + (exp(i?,(^V))/(a;)U=o. (5) 

where V — d/dx is one dimensional differential operator 
and the function f{x) is given by 

/(x) =tanh[/3(x-t-/i(i))], (6) 

as for that the effect of the differential operator on a 
function f{x), 

exp{a\7)f{x)\x=Q = f{x + a)\x=o, (7) 
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with any real constant a. By using the energy expres- where b is any real constant in a representative manner, 

sion given in Eq. Q , in Eq. ([5]) then, in order to get we write the exponential term in terms of the hyperbohc 

a polynomial form of the second term on the right-hand trigonometric functions, Eq. ([5]) exactly written in terms 

side of the equations, by using the van der Waerden iden- of muhiple spin correlation functions occurring on the 

tity for two-state spin, i.e., exp(6si) = cosh(6)-f sinh(6) right-hand side. Thus we get 

I 



,(1) 



)+{n[^i+4'^^i][^2+4''52]), 



\S=1 



(5=1 



(fe)\ 



,(fe-i) 



S3]n[^3 + .fB3][v43 + .sr^'i?3] , 



5=1 



(5=1 



(8) 



r 



with 



a decoupling approximation (DA) jSa] as 



An = cosh(J„V)/(a;)|a;=o, 
Bn = sinh(J„V)/(a;)|a;=o, 



where n = 1, 2, 3. 



is 



(k) Jfe) 



,(fe) 



In fact, the primitive form of this improvement corre- 
(^9) sponds essentially to the Zernike approximation [s^ in 
the bulk problem, and has been successfully applied to a 
great number of magnetic systems including the surface 
problems [40l - |43 |. By using the expansion in Eq. (fTO)) 
and the assumption below 



^nik, k = l,...,L, 



When the product in Eq. ([8]) is expanded, the multi 
site spin correlations appear. In order to make the ex- 
pansion manageable let us handle these correlations with the Eq. ([5]) morph into the form of, 

I 



(11) 



dm I 

dt 
dm2 

dt 



-mi -I- [Ai + miBiY[A2 -t- TO2S2], 

-1712 + [A2 + miB2][A3 + m2B^Y[A:i + rngSg], 



^cfTOfe ^ ^ [Az + mk-iB^WA'i + wikB^iY'lA'i + wik+iBz], 
dt 



drriL-i 
Jt 

drriL 

T ; 

dt 



-rriL-i + [A3 + mL-2B3][A3 + mL-iSg]^ A2 + miBa], 
-niL + [A2 + mL-iB2][Ai + ttilBiY ■ 



(12) 
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By using the binomial expansion and writing the hyper- bohc trigonometric functions in terms of the exponential 

functions we get the most compact form of Eq. ([T^ as 

I 



mi 



m2 



rrik 



niL = - 



z 1 \ 
-mi + ^ ^ Ai(7, rii)m1ml J , 

7=0t;=0 / 
z 1 1 \ 

-TO2 + ^ ^ ^ As (7, r/, i^)mlm1m'^ j , 

7=0 r;=0 1^=0 / 



-TOfc 



7=0r;=0i/=0 / 

z 1 1 \ 

-rriL-i + Z Z Z ^2(7, v)mlrnl_^m'ij^^ \ , 

7=0 i)=0 i/=0 / 
7=0>7=0 / 



(13) 



where 

Ai(7,?]) 
A2(7,?/,'^) 
A3(7,?/,'^) 



Ar-'Al-^BjBl 



(14) 



Note that, throughout our calculations t — 1 for simplic- 
ity. Self-consistent equations in Eq. (IT3| are typical first 
order ODE but has not an analytical solution because 
of the right-hand side contains transcendental functions. 
The dynamical equation of motion can be solved by var- 
ious numerical methods. In this work, we prefer to use 
the fourth order Runge-Kutta method (RK4) to get the 
evolution of the mit) by regarding Eq. (fT3|) as an ini- 
tial value problem. We can mention that the differential 
equation derived in Eq. (fT3)) extends up to the term 
. Each term in the equation of motion makes contri- 
bution to the solution, because the value of m, which is 
calculated at each time step, is iteratively related to the 
previous m value, however the situation is different from 
the behavior of the equilibrium systems at which high 
ordered terms can be neglected in the neighborhood of 
phase transition point. 

The system has three dependent Hamiltonian vari- 
ables, namely frequency of external magnetic field w, am- 
plitude and the film thickness L. For certain values 
of these parameters, temperature and the Ji, J2, J3 in- 
teraction constants, RK4 will give convergency behav- 
ior after some iterations i.e. the solutions have property 



I 

m{t) — m[t + 27r/w) for arbitrary initial value for the 
magnetization (m(i = 0)). Each iteration, i.e. the calcu- 
lation magnetization for t-\-\ from the previous magneti- 
zation for i, is now performed for these purpose whereby 
the RK4 iterative equation is being utilized to determine 
the magnetization for every i. In order to keep the itera- 
tion procedure stable in our simulations, we have chosen 
lO'* points for each RK4 step. Thus, after obtaining the 
convergent region and some transient steps (which de- 
pends on Hamiltonian parameters and the temperature) 
the layer's average magnetization can be calculated from 



Qk 



2ti 



■mk{t)dt 



(15) 



where is a stable and periodic function anymore and 
finally the DOP can be calculated by the arithmetic mean 
as 



1 ^ 



(16) 



k=l 



As cut-off condition for numerical self-consistency, we de- 
fined a tolerance 



< 10' 



(17) 



meaning that the maximum error as difference between 
the each consecutive iteration should be lower than 10~^ 
for all step. On the other hand, the modified surface 
exchange interaction has been defined to determine the 
different characteristic behavior of the system in certain 
range as 



Ji = Ml + A,). 



(18) 
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There are three possible states for the system, namely 
F, P and the coexistence phase (F+P). The total magne- 
tization time-series ■m{t) in convergent region is satisfied 
by this condition 

m{t) = -m{t + 7r/uj) (19) 

in the P phase which is called the symmetric solution. 
The solution corresponding to P phase follows the exter- 
nal magnetic field and oscillates around zero value which 
means that the DOP is zero. In the F phase, the solution 
does not satisfy Eq. ([TO]) and this solution is called as 
non-symmetric solution which oscillates around a non- 
zero magnetization value, and does not follow the exter- 
nal magnetic field i.e. the value of Q is different from 
zero. In these two cases, the observed behavior of mag- 
netization is regardless of the choice of initial value of 
magnetization m(0) whereas the last phase has magneti- 
zation solutions symmetric or non-symmetric depending 
on the choice of the initial value of magnetization corre- 
sponding to the coexistence region where F and P phases 
overlap. The main aim of the treatment is manifesting 
the frequency dispersion of the critical temperature co- 
ordinates of special point can be calculated by benefiting 
from the phase diagrams in (fc^T/ J3 — A^) planes for dif- 
ferent field amplitudes in order to understand and clarify 
the behavior of the dynamical system. 

III. RESULTS AND DISCUSSION 

Throughout our review, the best appropriate values 
and/or intervals of parameters are chosen and surface 
exchange Ji is directly modified on Ag by passing the 
distinction between bulk and interface exchange as J2 = 
J3 = 1.0. Unlike the deductive analysis, presenting the 
inferences as a result of investigation by induction is an 
essential principle in terms of understanding the overall 
behavior. In this sense, variation of DOP with tempera- 
ture for the selected films with different L = 4, 7 layers 
which has been rigorously depicted in the text, consti- 
tutes as the starting point for our systematic investi- 
gation, is presented in Fig. ([2]) for two selected values 
of external field amplitude as ho/Js = 1.0, 2.5 respec- 
tively. Also can be seen for the principal three frequency 
agents uj — 0.5, 1.0, 4.0. Incidentally, it is worth noting 
that the value of the crossover point A* has been calcu- 
lated as 0.3053 within the EFT as reported by Kaneyoshi 
[27j . A little bit the left and the right of this point at 
As = 0.0053, and 0.6053 respectively, corresponding av- 
erage layer magnetization has contra-arranged strength 
of the magnetic order mutually. Thus, the three repre- 
sentatives value of the modified exchange also has been 
taken and denoted from (a) to (b) as competence param- 
eter A, = 0.0053, 0.3053, and 0.6053. We also point 
out that this labeling procedure can be used for usual 
relevant behaviors in the same investigation from a dif- 
ferent viewpoint, i.e. it has been used for temperature 
dependency of the average magnetization of each layer 
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FIG. 2. (Color online) Variation of dynamical order pa- 
rameter with temperature for the selected films with differ- 
ent L = 4, 7 layers for two values of external field ampli- 
tude. The left-hand- and right-hand-triple-panels correspond 
to ho/Js = 1.0 and 2.5 respectively can be seen for the prin- 
cipal three frequency agents as a; = 0.5, 1.0 and 4.0 stated 
in the text. The number accompanying each set of curves 
denote the values of modified exchange from (a) to (b) as 
A, = 0.0053, 0.3053, 0.6053. 



in Fig. ^ and the variation of it with the layer index 
/c at a fixed temperature where ksT/ J3 — 3.5 in Fig. 
(|4]) . Average magnetization profiles across the film differ 
qualitatively in two regimes (As < A* and A^ > A*) as 
shown in Fig. ^ for the amplitude agents. Apart from 
this, DOP and consequently the critical temperature is 
independent from thickness L in the crossover point and 
this makes it 'special'. The thickness-independent crit- 
ical temperature of the film at special point reduces to 
the bulk with coordination number z = 6 critical value 
ksT/ Jt, — 5.0734 in infinite-frequency for particular am- 
plitude value or zero-amplitude limit while there is no 
significant change in the value of the A* when w or ho 
chances. This reduction is one of the major indicators 
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FIG. 3. (Color online) Variation of dynamical order parame- 
ter for each layers with temperature for the selected films with 
different L — 4, 7 layers for two values of external field ampli- 
tude. The interchange of the layer hierarchy can be monitored 
in the figure for the same values depicted in Fig ((2]). 
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FIG. 4. (Color online) The average magnetization profiles for 
a film L = 7 layers for two values of external field amplitude 
ho/Js — 1.0, 2.5 and fixed temperature ksT/Jz = 3.5. The 
number accompanying each curve denotes several values of 
modified exchange from (a) As = 0.0053 to (b) As = 0.6053. 
The evolution of the curves can be seen in the left-hand- and 
right-hand-triple-panels for the values of frequency w = 0.5, 
1.0 and 4.0. 



for the accuracy of our treatment. The effect of ampli- 
tude or frequency of the external oscillatory magnetic 
field on critical behavior of the bulk systems is now well 
known. At first sight, we can see in Fig. ([2]), ([3]), and 
Fig. (HI that a weak increasing of w causes an increase 
in fcsTc/Js. However, an increasing of hg/Js causes a 
decrease ksTc/Js for given and arbitrary uj. This 
is an expected result, since the increase in frequencies 
enhances the phase lag between time dependent mag- 
netization and the external field signals. Therefore, the 
asymmetric behavior of HL in the form of Lissajous curve 
which is m(t) of the system versus the external magnetic 
field h{t), becomes more prominent, since the time depen- 
dent magnetization has less time to follow the oscillatory 
field. Hence, the system can undergoes a DPT which re- 
quires a small amount of thermal energy. On the other 
hand, according to times-series of the magnetization and 



the external magnetic field, increasing the field frequency 
at first, obstructs the saturation of the ordinary magne- 
tization due to the decreasing energy coming from the 
oscillating magnetic field in a half-time period which fa- 
cilitates the late stage domain growth by tending to align 
the moments in its direction (i.e. the magnetization be- 
gins to fail following the oscillatory field) and this makes 
the occurrence of the frequency increasing route to DPT 
at the critical point. There is a concurrence by different 
researchers that the HL loses its symmetry when the os- 
cillating period of external perturbation becomes much 
smaller than the typical relaxation time of the system. 

The above-mentioned mechanism is due to the com- 
peting time scales in such non-equilibrium systems and 
also can be seen in the global phase diagrams in {ho/ — 
ksT/ J3) planes presented in Fig. ([5]). The external oscil- 
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latory field defined in Eq. ([2]) is governed by only an open 
variable /ig, which is the amplitude of the field. However, 
the system has close dependency to frequency since it is 
embedded into a transcendental function. Although the 
time average of the field over a full period is zero, sinu- 
soidally oscillates field drag the system to the ordered 
phase. In addition to this, bulk and modified exchange 
interactions enforce the system to stay in the ordered 
phase. The last factor, the temperature, which causes 
thermal agitations, induces a disordered phase when the 
energy supplied by the temperature to the system is high 
enough. Thus, mainly a competition takes place between 
the amplitude and temperature. After a brief summary 
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FIG. 5. (Color online) The phase diagram in the {ho/Ja — 
ksT/Jz) plane for the films with difi^erent L = 4, 7 lay- 
ers. Frequency changes with uj = 0.1,0.25,0.5,1.0,2.0 and 
4.0 for each panels. The numbers accompanying each set of 
curves denote the value of modified exchange from (a) to (b) 
as As = 0.0053, 0.3053, 0.6053 and this type of display format 
is also stated as an arrow in the figure. The dashed lines in 
the middle-left-hand panel corresponds to the representatives 
serve as a viewing guide. 



of the physical background of the critical phenomena ob- 
served in the dynamical system, let us investigate the 
effects of external field frequency and modified exchange 
on the selected films of different L = 4, 7 thickness. In 
Fig. ([5]), phase diagrams can be seen for the relevant 
L with chosen A^. At first glance, one can see that 
the Ising-type pure crystalline thin film exhibits meta- 
stable behavior such as storable coexistence region at low 
temperature and high amplitude values where the F and 
P phases overlap and the critical properties of the sys- 
tem depends on the initial value of the magnetization. 
If the frequency value increases such meta-stable phases 
appears. Contrary to this, the coexistence (F-l-P) regions 
disappear with decreasing frequency at the high temper- 
ature and low amplitude values. Although there is a 
noteworthy information about the meta-stability, proba- 
bly it would be wrong to skip the unpretentious remark 
that the meta-stable regions also have the maximal area 
in crossover value of the modified exchange parameter in 
the diagrams. The representative amplitudes has already 
been selected as h^/Js, = 1.0 and 2.5 privately to avoid 
complications of coexistence regions that may be an arti- 
fact of the effective-field approach itself reported several 
times by some researchers in such dynamical phenomena 
[43l - l46j . In addition, there exists a dynamical tricriti- 
cal point on the dynamical phase boundary at which the 
second order and first order phase transition lines are 
separated from each other. Besides, rising frequency led 
to the flaring curves in the phase planes, contrary to de- 
creasing frequency shrinks the ferromagnetic region. This 
reflects the truth, since increasing the field frequency 
causes a growing phase delay between the magnetization 
and held, then this makes the occurrence of the DPT dif- 
ficult, as a result of this mechanism the dynamical phase 
boundary (DPB) gets wider. The modified exchange has 
been chosen as the best appropriate values A.; = 0.0053, 
0.3053, 0.6053 with reference to above procedure (which 
provide the conditions A^ < A* and A* < Ag, respec- 
tively) and it is clear that the results presented in Fig. 
([5|) are completely consistent with the results shown in 
Fig. (21), (121), and Fig. (UJ since the Ag comes close to 
the crossover value then the behavior of the system be- 
gins to resemble as quasistatic (the phrase is coming from 
the fact that the HLs has a nonzero residual loop area, 
originates from the deficiency of the relevant method, 
even in the zero- frequency limit). Also, the interchang- 
ing phenomenon can be seen from the critical value of the 
transition in inner- and outer-stablemate-curves around 
the characteristic one corresponding crossover. Physical 
mechanism mentioned above can be briefly explained as 
follows: If one keeps the system in one well of a Landau 
type double well potential, a certain amount of energy 
originating from magnetic field is necessary to achieve 
a dynamic symmetry breaking. If the amplitude of the 
applied field is less than the required amount then the 
system oscillates in one well. In this situation, the mag- 
netization does not change its sign. In other words, the 
system oscillates around a nonzero value. This region is 
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dynamically ordered phase. When the temperature in- 
creases, the height of the barrier between the two wehs 
decreases. As a resuh of this, the less amount of mag- 
netic field is necessary to push the system from one well 
to another and hence the magnetization can change its 
sign for this amount of field. Consequently, the time 
averaged magnetization over a full cycle of the oscillat- 
ing field becomes zero. In addition, we can mention 
that for a dynamical thin film, increasing the dipole- 
dipole interaction-induced energy contribution by rising 
the strength of the modified surface exchange interaction 
will cause an increase the orderliness tendency of the sys- 
tem which makes the transition from one well to another 
more difficult. 

As stated in various place, a prime objective of the 
study was to confirm whether a common value of ex- 
isted for all films, and if so to determine its value. This 
requires a very accurate determination of the transition 
temperature for each film thickness and each value of sur- 
face exchange enhancement. For this purpose, we need to 
dart a glance at the relevant (Ag — ksT/ J3) global phase 
diagrams, thus obtained, are shown in Fig. (|6]) which are 
plotted for different numbers of layers when the strength 
of the field amplitude is h^/Jz — 1.0, and 2.5 respec- 
tively. As As rises, the critical temperature of the system 
increases for arbitrary L. For each value of the frequency 
with the representative amplitudes, as seen clearly that 
the curves with different film thicknesses intersect each 
other at the same abscissa point A^ = A* = 0.3053. 
According to this result, the fact that the special point 
can be defined as that particular Aj value at which bulk 
critical temperature is independent of the film thickness 
L occurs, it is assumed to coincide with the critical tem- 
perature ksT^ /J3 of the usual simple cubic lattice. Fur- 
thermore, based on the definition of A^, it can be ex- 
pected that the crossover point in Fig. ([6]) should define 
also the critical temperature of three-dimensional infinite 
bulk system, where the surface and the A^ parameter are 
of no importance. This is really the case, which can be 
seen also from Fig. ([B]), where the bulk ksT^ / and the 
surface ksT^ / J3 critical temperatures of the correspond- 
ing semi-infinite system are represented, respectively, by 
the dashed lines and also pointed by one each arrow. For 
the values of Ag < A*, thin films have lower critical 
temperature than the bulk system. The value of A^ also 
changes the relation between the thickness and the criti- 
cal temperature of the film. ksTd increases with the 
film thickness L, and approaches ksT^ / Jz asymptoti- 
cally as the number of layers become large. The thicker 
films have higher critical values for As < A* whereas 
they have lower critical values for A^ > A* than thinner 
ones. It is worth noting once again that the critical tem- 
perature of the film is independent of L when A^ — A* , 
and equal to ksT^ I Jz- On the other hand, for Aj > A*, 
the critical temperature ksTc/ J3 is greater both than the 
bulk ksT^ I J3 and the surface ksT^ /J3 critical tempera- 
tures of the corresponding semi-infinite Ising system and 
larger the L is, the lower ksTc/Js is. The film criti- 
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FIG. 6. (Color online) The critical temperature versus mod- 
ified exchange in (As — fcsT/Ja) phase planes for the films 
with different L = 4, 5, 6, and 7 layers. Two set of curves 
correspond to the selected amplitude agents /10/J3 = 1.0 and 
2.5. The number accompanying each panels denote the value 
of frequency uj — 0.1,0.25,0.5,1.0,2.0 and 4.0. Differences 
in the layer hierarchy on both sides of the crossover point is 
shown by the arrows in opposite directions. The exemplary 
inset of the crossover point is stated in terms of visual clarity. 



cal temperature A:bTc/ J3 approaches asymptotically the 
surface critical temperature fc^T^/ J3 of the correspond- 
ing semi-infinite system as the number of layers become 
large. Beyond all results, the evolution of the curves with 
respect to frequency as w = 0.1, 0.25, 0.5, 1.0, 2.0, and 
4.0 for particular two amplitude agents as stated several 
times above IiqI = 1.0, and 2.5 can observed in Fig. 

Accordingly, causes all the curves go- up and -down 
with rising ui/hQ/J^ respectively also the same opera- 
tion increases/decreases the critical temperature of the 
system while it does no affect the relation between the 
film thickness and critical temperature, i.e. for Ag < A* 
thicker films have higher critical values and reverse is 
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FIG. 7. (Color online) The frequency dispersion of the critical 
temperature coordinate of the special point for different values 
of the amplitude of external oscillatory field ho/J-j, — 0.0, 
0.25, 0.5, 1.0 and 2.5. The dashed line at cj = 0.1 serves as a 
boundary for quasistatic limit. 



valid for As > A* with ui, ho ^ 0, if the special point 
present. In addition, there is no significant change in the 
abscissa point of the crossover (A*) while rising or lower- 
ing of two oscillatory parameter transports the phase di- 
agrams collectively in the space. When lu is large enough 
for finite amplitude value as ho ^ 0, the ordinate value 
of the special point approaches the bulk critical temper- 
ature in quasistatic case ksT^ /J^ asymptotically. The 
observations about dynamical feature of the global phase 
diagrams in (Ag — ksT/J^) planes can also be explained 
by the well known mechanism underlying the dynamical 
phase transitions phenomena. In other words, the afore- 
mentioned scenario is onset as things stand. 

As one of the most striking aspects of our work is that 
the frequency dispersions of the critical temperature co- 
ordinates of special point for various external oscillatory 
field amplitudes as ho/Js = 0.0, 0.25, 0.5, 1.0, 2.5 shown 
in Fig. ([7]). In accordance with our anticipations, the dis- 
persion curves evolve out of the characteristic line which 
accounts for the pure bulk system in quasitatic limit as 
rising ho/J^. The behavior of the dispersions qualita- 
tively supports the last parts of the descriptions in the 
previous paragraph that is to say all dispersion curves 
saturates to the value of bulk critical temperature in 
w — >■ oo limit. Contrary to this, they goes to a finite 
non-zero value in u ^ limit. Major attention has been 
paid to catch the overall behavior of a dynamical thin 
film in our calculations, hence the lowest value of the 
frequency in the treatment is lo = 0.1. This value is suf- 
ficient for estimation an existence of the HL residuality 
mentioned above even now. 



IV. CONCLUSION 

In this study, we have investigated the surface enhance- 
ment phenomenon and dynamic nature of the critical 
phenomena which is observed for dynamical Ising-type 
thin films with various thickness which has been defined 
with inner coordination number z = 4 driven by an exter- 
nal oscillatory magnetic field by means of effective-field 
theory based on a standard decoupling approximation 
and the time evolution of the system has been presented 
by utilizing a Glauber type stochastic process. Our start- 
ing point for the systematic review was to examine the 
DOP, average magnetization of layers with respect to 
temperature and layer index for the selected field am- 
plitude and frequency values at finite temperature below 
the Curie point of the static case on a pure Ising-type 
bulk. In addition to the observations set forth in the 
literature, the effect of the frequency and amplitude on 
the standard arguments have been propounded. For this 
purpose, the best appropriate values and/or intervals of 
parameters which are given in the text, have been chosen 
for the treatment. We supported the relevant investiga- 
tions with phase diagrams globally as well as the average 
magnetization. 

In order to examine the effect of the dynamical param- 
eters of external field on the phase diagrams, firstly we 
presented the phase diagrams for selected three distinct 
modified exchange representatives with three values of 
frequency in (/lo/ J3 — ksT / J'^) planes. In the immediate 
aftermath, the frequency induced evolution of the cor- 
responding agents along the routes at ho/J^ — 1.0, 2.5 
respectively in the previous phase diagrams have been 
investigated at finite temperature in (Ag — ksT/Js). 
According to our findings, a weak decreasing of the 
amplitude causes decreasing in A^ coordinate of the 
crossover point, while a weak increasing of the frequency 
causes an increasing critical temperature component of 
the crossover. For a better view, we have depicted fre- 
quency dispersion of the ksTc/ J3 coordinate of the spe- 
cial point for aforementioned amplitude agents. 

EFT takes the standard MFA predictions one step for- 
ward by taking into account the single spin correlations 
which means that the thermal fiuctuations are partially 
considered. Although all of the observations reported 
in this work shows that EFT can be successfully ap- 
plied to such nonequilibrium realistic systems, the true 
nature of the physical facts underlying the observations 
displayed in the system (especially the origin of the differ- 
ent contrary effects of two variables of the external filed 
related between the competing time scales) may be fur- 
ther understood with an improved version of the present 
EFT formalism which can be achieved by attempting 
to consider the multi spin correlations which originate 
when expanding the spin identities or a MC simulation 
as stochastic process. We believe that this attempt could 
provide a treatment beyond the present approximation. 

In conclusion, we hope that the results obtained in 
this work would shed light on the further investigations 
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of the dynamic nature of the critical phenomena in pure 
crystaUine systems (e.g. thin films and semi-infinite sys- 
tems) and would be beneficial from both theoretical and 
experimental points of view. 

ACKNOWLEDGEMENTS 

The numerical calculations reported in this paper were 
performed at TUBiTAK ULAKBIM (Turkish agency), 



High Performance and Grid Computing Center (TRUBA 
Resources) and this study has been completed at Dokuz 
Eyliil University, Graduate School of Natural and Ap- 
plied Sciences. One of the authors (B.O.A.) would like 
to thank the Turkish Educational Foundation (TEV) for 
full scholarship. 



D. Weller, S. F. Alvarado, W. Gudat, K. Schroder, and 
M. Campagna, Phys. Rev. Lett. 54, 1555 (1985). 

W. Diirr, M. Taborelli, O. Paul, R. Gcrmar, W. Gu- 
dat, D. Pescia, and M. Landolt, Phys. Rov. Lett. 62, 206 
(1989). 

Yi Li and K. Babcrschkc, Phys. Rev. Lett. 68, 1208 
(1992). 

M. Farle, K. Babcrschkc, U. Stetter, A. Aspelmeier, and 
F. Gerhardtcr, Phys. Rev. B. 47, 11571 (1993). 
P. Poulopoulos and K. Baberschke, J. Phys.: Condens. 
Matter. 11, 9495 (1999). 

B. T. Jonker, K. H. Walker, E. Kiskcr, G. A. Prinz, and 

C. Carbone, Phys. Rev. Lett. 57, 142 (1986). 

D. Pescia, M. Stampamoni, G. L. Bona, A. Vaterlaus, R. 
F. Willis, and F. Meier, Phys. Rev. Lett. 58, 2126 (1987). 
D. T. Pierce, Surf. Sci. 710, 189 (1987). 

D. P. Pappas, K. P. Kamper, and H. Hopster, Phys. Rev. 
Lett. 64, 3179 (1990). 

C. Ran and S. Eichncr, Phys. Rev. Lett. 47, 939 (1981). 
C. Rau, C. Jin, and M. Robert, J. Appl. Phys. 63, 3667 
(1988). 

M. Stampanoni, A. Vaterlaus, M. Acschlimann, and F. 
Meier, Phys. Rev. Lett. 59, 2483 (1987). 
M. Przybylski and U. Gradmann, Phys. Rev. Lett. 59, 
1152 (1987). 

C. M. Schneider, P. Bressler, P. Schuster, J. Kirschner, 
J. J. de Miguel, and R. Miranda, Phys. Rev. Lett. 64, 
1059 (1990). 

Th. Dctzcl, M. Vonbank, M. Donath, N. Mcmmel, and 

V. Dose, J. Magn. Magn. Mater. 152, 287 (1996). 

M. Stampanoni, Appl. Phys. A. 49, 449 (1989). 

A. S. Edelstein, C. Kim, S. B. Qadri, K. H. Kim, V. 

Browning, H. Y. Yu, B. Maruyama, and R. K. Everett 

Solid State Commun. 76 1379, (1990). 

R. Krishnan, H. O. Gupta, H. Lassri, C. Sella, and M. 

Kaabouchi, J. Appl. Phys. 70, 6421 (1991). 

E. Colombo, O. Donzelli, G. B. Fratucello, and F. Ron- 
coni, J. Magn. Magn. Mater. 107, 1857 (1992). 

Yi Li, C. Polaczyk, F. Kloso, J. Kapoor, H. Malctta, F. 
Mezei, and D. Ricgel, Phys. Rev. B 53, 5541 (1996). 



[21 

[22; 

[23; 

[24' 
[25^ 

[26 

[27 

[28; 

[29; 

[30 

[31 

[32 

[33 
[34 
[35^ 
[36^ 

[37 
[38^ 
[39^ 
[40^ 
[41 

[42; 

[43; 

[44 

[45; 

[46; 



M. Farle and K. Baberschke, Phys. Rev. Lett. 58, 511 
(1987). 

D. Pescia, G. Zampieri, M. Stampanoni, G. L. Bona, R. 
F. Willis, and F. Meier, Phys. Rev. Lett. 58, 933 (1987). 
K. Binder and P. C. Hohenberg, Phys. Rev. B 9, 2194 
(1974). 

K. Binder, Thin Solid Films 20, 367 (1974). 

F. Aguilera-Granja and J. L. Moran-Lopez, Phys. Rev. 

B 31, 7146 (1985). 

D. P. Landau and K. Binder, Phys. Rev. B 41, 4633 

(1990) . 

T. Kaneyoshi, J. Phys.: Condensed Matter 3, 4497 

(1991) . 

Y. G. Wang, W. L. Zhong and P. L. Zhang, SoUd State 
Commun. 101, 807 (1997). 

D-L Yao, Y-Z Wu, and Z-Y Li Phys. Stat. Sol. B, 231, 
3 (2002). 

M. Saber, A. Ainanc, F. Dujardin, and B. Stebe, Phys. 
Rev. B 59, 6908 (1999). 

Y. Laosiritaworn, J. Poulter, and J. B. Staunton, Phys. 
Rev. B 70, 104413 (2004). 

P. Cossio, J. Mazo-Zuluaga, and J. Restrepo, Physica B 
384, 227 (2006). 

Y. Laosiritaworn, Thin Solid Fihns 517, 5189 (2009). 
U. Akmci, J. Magn. Magn. Mater. 329, 178 (2013). 
R. J. Glauber, J. Math. Phys. 4, 294 (1963). 
R. Honmura and T. Kaneyoshi, J. Phys. C: Solid State 
Physics 12, 3979 (1979). 

T. Kaneyoshi, Acta Phys. Pol. A 83, 703 (1993). 

J. W. Tucker, J. Magn. Magn. Mater. 102, 144 (1991). 

F. Zernike, Physica A 7, 565 (1940). 

T. Balcerzak, J. Magn. Magn. Mater. 97, 152 (1991). 
T. Kaneyoshi, R. Honmura, I. Tamura, and E. F. Sar- 
mento, Phys. Rev. B 29, 5121 (1984). 
T. Kaneyoshi, I. Tamura, and E. F. Sarmento, Phys. Rev. 
B 28, 6491 (1983). 

B. O. Akta§, U. Akmci, and H. Polat, Physica B 407, 
4721 (2012). 

O. Idigoras, P. Vavassori, and A. Berger, Physica B 407, 
9 (2012). 

G. Bcrkolaiko and M. Grinfeld, Phys. Rev. E 76, 061110 
(2007). 

X. Shi and G. Wei, Phys. Lett. A 374, 1885 (2010). 



